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In General Relativity, the constraint equation relating metric and density perturbations is inherently nonlinear, 
leading to an effective non-Gaussianity in the dark matter density field on large scales - even if the primordial 
metric perturbation is Gaussian. Intrinsic non-Gaussianity in the large-scale dark matter overdensity in GR is 
real and physical. However, the variance smoothed on a local physical scale is not correlated with the large- 
scale curvature perturbation, so that there is no relativistic signature in the galaxy bias when using the simplest 
model of bias. It is an open question whether the observable mass proxies such as luminosity or weak lensing 
correspond directly to the physical mass in the simple halo bias model. If not, there may be observables that 
encode this relativistic signature. 


I. INTRODUCTION 

In Newtonian gravity, the Poisson equation is a linear relation between the gravitational potential and the matter overdensity. 
By contrast, in General Relativity (GR) this is replaced by a nonlinear relation, which introduces mode coupling between large 
and small scales lUHSl. The original result has been confirmed by a number of independent calculations mil. Similar mode 
coupling can be produced in Newtonian gravity by local-type primordial non-Gaussianity of the gravitational potential mniin]. 
The GR effect has therefore previously been interpreted as an effective local non-Gaussianity on very large scales HEIEHia. 

Recently two papers have argued that a “separate universe” approach can be used to show that no scale-dependent bias 
arises from the GR corrections on large scales cmni. The argument is that the nonlinear coupling between long-wavelength 
perturbations on a scale Xp, and the small-scale variance, cr| = (<5|), on a scale Ag, vanishes under a local coordinate rescaling 
and hence is unobservable. 

The separate universe approach GSlIIll has proved to be a powerful tool to understand the origin of large-scale structure, 
and primordial non-Gaussianity, from inflation. Accelerated expansion in the very early Universe stretches initial small-scale 
vacuum fluctuations up to scales much larger than the Hubble scale at the end of inflation. Spatial gradients for such long- 
wavelength modes become small relative to the local Hubble time, and for many scales of interest, the perturbed universe can 
be treated as a patchwork of “separate universes”, each locally obeying the classical Friedmann-Lemaitre-Robertson-Walker 
(FLRW) evolution of an unperturbed universe. 

The separate universe approach is particularly useful for studying nonlinear perturbations on large scales |[T8l|20|. For adia¬ 
batic perturbations, each separate universe patch follows locally the same evolution as the unperturbed “background” cosmology. 
The only difference between separate patches is the local expansion, characterised by the comoving metric perturbation C,. This 
is defined to be the local perturbation of the integrated expansion rate with respect to a background flat reference cosmology, 
5N = N — N, where N = J dtQ/3. 

An important consequence of the uniqueness of the local evolution for adiabatic perturbations is that ( is conserved on large 
scales where the separate universe approach is valid mmiiiaii. 
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FIG. 1: Schematic of the various scales in 0. 


In each patch, the comoving spatial line element is (see ini) 

ds(3) = Sijdx^dx^ . (1) 

There is a global background which must be dehned with respect to some scale Aq, at least as large as all the other scales of 
interest, i.e., at least as large as our presently observable Universe. It is important to distinguish this from the scale of the separate 
universe patches, Ap. This is large enough for each patch to be treated as locally homogeneous and isotropic, but patches must 
be stitched together to describe the long-wavelength perturbations on a scale Ap ^ Ap. Thus, following in, we require a 
hierarchy of scales (see Fig. 1); 


Aq > Ap ^ Ap ^ As . (2) 

The local observer in a separate universe patch cannot observe the effect of ^p, which is locally homogeneous on the patch 
scale Ap. However, local coordinates can be dehned only locally and the long mode curvature perturbation is observable through 
a mapping from local to global coordinates. 


II. THE PHYSICAL EFFECT OF CURVATURE WITHIN THE OBSERVABLE UNIVERSE 

In Newtonian gravity the only constraint on initial conditions is the Poisson equation, which provides a linear relation between 
the overdensity and the gravitational potential at all orders 

= (3) 

Thus if the initial Newtonian potential is Gaussian, then so is the initial density held 6. In GR, the nonlinear energy constraint 
equation for irrotational dust is f23\ 


-2a^ + = 16 ttGp + 2A, (4) 

o 

where p is the comoving matter density, A is the cosmological constant, 0 = is the expansion rate of the matter 4-velocity, 

a is its shear, and is the Ricci curvature scalar of the 3-dimensional space orthogonal to it^. At hrst order in perturbations 
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about an FLRW cosmology, the energy constraint combines with the momentum constraint to give the relativistic version of 
the Poisson equation ([^, where is replaced by $, i.e. the spatial metric perturbation in longitudinal gauge, and 5 is the 
synchronous comoving gauge density contrast. Note that $ = 3^/5 at first order. At second order, at the start of the matter era, 
using the relation between and (^, we obtain m 

= (5) 

Consider a Gaussian distribution of (. We separate C and 6 into independent long- and short-wavelength modes, C = Cl + Cs 
and S = Sl + Ss, where the wavelength of the long modes Xl obeys in particular, Xl ^ Xp. To leading order in Cs and Cl, 
and neglecting gradients of Cl relative to those of Cs, the initial constraint (m implies V^Cl = and 

T/\s-2(:L^\s = -la^HHs- ( 6 ) 

The second term on the left represents the long-short mode coupling. 

Within a local patch on a scale Xp, it is possible to redefine the background spatial coordinates to absorb the effects of the 
long-wavelength perturbations Cl, following iflTl : 

i^ = x^+e, e = CLx\ (7) 

If we neglect gradients of the long mode, this transformation eliminates Cl from the spatial metric Q 

'^^( 3 ) ~ 5ijdx'‘dx^ = Sijdx^dxR ( 8 ) 

This transformation holds for each single patch (see Fig. 1). 

Since this is a purely spatial coordinate transformation, the curvature and density perturbations transform as scalars, 

(six) = Cs{x ), h{x) = Ss{x), where x = [l + C,l{x)\x . (9) 

The constraint equation (j^ becomes 

V^Cs{x) = -\a^HHs{x). (10) 

Thus in the new local coordinates, in one patch of size Xp, we have a linear Poisson equation and the long-short mode coupling 
appears to be absent. This confirms the fact that the local observer in a separate universe patch cannot observe the effect of the 
locally homogeneous perturbation Cl, as argued in ifT^fTTll . 

The original coordinates x* define a global chart, which is essential for defining random fields such as Cl on large scales, 
and the long mode curvature perturbation enters through the mapping from local to global coordinates. Indeed, a coordinate 
transformation that depends on a random field is not a new concept in large-scale structure. The situation here is reminiscent 
of the redshift-space distortion map, where the random field is given by the peculiar velocities (which are in turn generated by 
large-scale density perturbations). Because of the nonlinear nature of this map, an initially Gaussian field in real space becomes 
non-Gaussian in redshift space Il24ll25l . 

As shown in iflTl . C(i) and 5{x) are Gaussian fields with respect to the local coordinates i®. By (|^, 5{x) = 5{x) and so we 
recover a non-Gaussian distribution for the density field with respect to the global coordinates, x®. See Fig. 2 for a schematic 
illustration of this. 

There is no argument against the fact that a;® coordinates are useful to discuss physics in a local patch of size Ap <C Al - but 
this applies only locally. The effect of the long mode C,l is to create spatial curvature V^Cl on a constant-time hypersurface. 
The curvature can be eliminated only locally, by neglecting gradients of C,l, as in ([^. Beyond the single patch, when gradients 
are not negligible, we have 1261 


= e^^^e^^^Sijdx^dx^ = 


Sijdx^dx^ + 0{\xfV^CL) 


( 11 ) 


Indeed, all coordinate transformations that neglect curvature can only be defined locally (see 1261 ). The spatial curvature 
generated by (p is directly related to the long-wavelength density perturbation Sp by the long-wavelength part of (|^. On a 
constant-time hypersurface, the density perturbation cannot be eliminated by the spatial coordinate transformation. This implies 
that we need many different local patches of scale ~ Xp described by different i®-coordinates in the entire observed universe of 
scale ~ Aq. 

Within a single patch on a scale Xp, the local observer does not notice this stochastic nature of the local coordinates. However, 
once we are interested in physics beyond the local patch we have to compare the different patches and we notice that i® vary 
stochastically through their dependence on This implies that GR non-Gaussianity is present in the dark matter distribution. 
Nonlinear GR effects can be appreciated only by looking at a region where V^(p is not negligible and, in principle, the stochastic 
nature of (p could become apparent. 
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FIG. 2: Illustrating the non-Gaussianity of Ss{x) and its correlation with in the global coordinates, x. We start from the fact that Ss{x) 
is Gaussian (HI, and uncorrelated with the Gaussian (,l{x). The top 2 panels show one k mode for each. Then we apply the transformation 
X ^ X = x[l -\- The resulting Sg(x) field (next panel down) is clearly modulated by (l- To see that it is non-Gaussian, consider 

another k mode, Ss 2 {x) (fourth panel). Clearly Ss 2 {x) is highly correlated to 5s{x), i.e., there are phase correlations. The local variances of 
these two modes at fixed coordinate scale are correlated with Cl (bottom panel). 
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III. LOCAL-TYPE NON-GAUSSIANITY IN THE NEWTONIAN DENSITY FIELD 

In Newtonian gravity, it is standard to parametrise primordial non-Gaussianity of local type in the Newtonian potential by 

= </5 + /nl - (<P^)) , (12) 

where yi is a Gaussian random field. If we split the Gaussian field into long and short modes, ip = ips ‘Pl, and use the same 
assumptions that lead to (|^, then the Poisson equation ([^ yields 

(l + 2/NLV5L)vVs = -^a2ij2jg. (13) 

This leads to a scale-dependent galaxy bias oniin]. 

Single-field, slow-roll inflation generates an almost Gaussian distribution for which remains Gaussian for adiabatic pertur¬ 
bations on super-Hubble scales through to the start of the matter-dominated era. For this Gaussian case, using the hrst-order 
relation C = 5(/?/3, the GR second-order constraint equation (j^ gives 

{l-^PL)v^ps = -la^H^Ss. (14) 

This has been used to argue that GR corrections to the comoving density held will also lead to a scale-dependent bias, with an 
effective (local) non-Gaussianity parameter EiinHia: 


rGR.eff 

/nl 


(15) 


While the Newtonian potential satishes the linear Poisson equation at all orders, in GR C obeys the nonlinear constraint (j^. 
We could thus dehne primordial non-Gaussianity of the comoving density held, e.g., from inhation, in terms of the primordial 
metric perturbation, Local-type primordial non-Gaussianity can be parametrised as 


c = 


P + /nl - {p^)) 


(16) 


Thus, the non-Gaussianity present in the comoving density held Q is similar to a Newtonian density held with local-type 
non-Gaussianity parameter 


.peff _ £ 

/nl ~ /NL 


5 

3 ■ 


(17) 


However, as we shall show in the next section, this Newtonian description of the density held in terms of an effective /nl is 
not sufficient to describe local halo abundance in GR. When considering halos, smoothing of the dark matter held must be 
introduced. The variance of 5{x) smoothed on a physical scale is not correlated with the large-scale modes Cl, unlike the case 
with primordial non-Gaussianity of local type. 


IV. SCALE-DEPENDENT BIAS AND SINGLE-FIELD INFLATION 

The presence of a large-scale perturbation can change the local abundance of collapsed dark matter halos either through a 
perturbation in the local “background” density, 5l, or through a perturbation in the local variance, (Jr ^ + SaR, on a 

comoving scale R ^ Ap, where we use a bar to denote quantities in the global background. 

For a Gaussian density held the small-scale variance is independent of the long-wavelength perturbations, but for a non- 
Gaussian held of the local-type in Newtonian gravity, as given in ( [T3] l, the small-scale variance is modulated by long-wavelength 
perturbations and we have 


4l^, = (l+4/NL<^i)4- (18) 

Here cr^ is the variance of the density held smoothed on the comoving length scale R, corresponding to a given physical mass 
in the global background. 


M = 


4 

3 


ttR^P . 


(19) 
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We have 5an cx /nlV^l oc f-NL5L/k‘j^ and hence local-type non-Gaussianity can lead to a strongly scale-dependent bias for 

halos nniini. 

In GR the long-wavelength metric perturbation, (l, also leads to a change in the variance at a fixed comoving radius, R. The 
density perturbations obey the nonlinear constraint equation. This nonlinearity modulates the density field by (l- 

= (1-2Cl)<5(x), (20) 

where S{x) refers to the density perturbations in the absence of while denotes that in the presence of (^l- From (j^ we 
have 

41^^ = (i-4a)4- (21) 

However in GR we must also account for the perturbation in the mass enclosed within a comoving radius R due to a long- 
wavelength metric perturbation 


M= ^7r(l + 3CL)i?"p. (22) 

To compare the local abundance of halos of a fixed physical mass, M in ®, in the presence of a long-wavelength perturbation, 
Ql, we must compare properties of the density field on a fixed local physical scale {1 + (l)R = R. Hence 


R = {1-Cl)R- (23) 

Note that we neglect the perturbation in the density, p, in ( [22| , since is negligible relative to for long-wave length pertur¬ 
bations. 

For a scale-invariant distribution of the curvature perturbation, C„ the variance of the density field at two different comoving 
scales Ri and i ?2 are related by 


= 



(24) 


for scales greater than the matter-radiation equality horizon scale, Ri, R 2 > 


Aeq. Using (23 1 , we obtain 


4 — (1 + 4Ci) 4 ■ 


(25) 


Thus combining ( |2T| ) and ( |25| ), we obtain 

‘^R\a=^R- ( 26 ) 

As shown in ( [T0] i, the small-scale density at a fixed local physical scale is independent of the long-wavelength perturbation. 
Thus the long-wavelength mode has no effect on the small-scale variance of the density field smoothed on a fixed mass scale. 

Although there is a change in the local variance at a fixed comoving coordinate radius, R, this is exactly compensated by the 
change in the local coordinate scale corresponding to a fixed physical mass and the scale-dependence of the variance ini. This 
is consistent with the argument that the long-wavelength perturbations Xl ^ Ap are not observable locally if the distribution of 
the primordial metric perturbation, C,, is a Gaussian random field imini. 

Note that in Newtonian gravity there is no equivalent perturbation in the physical volume and hence there is no perturbation 
in the mass enclosed at a coordinate radius R caused by the large-scale potential (/jp. A Newtonian description with = 

—5/3 thus fails to account for the local coordinate invariance inherent in GR. 


V. CONCLUSIONS 

We have shown how non-Gaussian correlations in the matter overdensity arise due to nonlinear constraints in GR, even when 
the primordial metric perturbation from inflation, (^, is described by a Gaussian random field. This may be understood in a simple 
way as arising from the long-wavelength metric perturbation rescaling the local small-scale curvature perturbation Cs, and 
thus the local small-scale density field, through the nonlinear constraint equation (j^. 

A very long wavelength metric perturbation can be removed locally by a coordinate transformation (0. If the long-wavelength 
perturbation were much larger than our observable horizon, Ap ^ Aq, then that would be the end of the story - perturbations 
with wavelength much larger than our horizon form part of our background cosmology and cannot be observed locally. However, 
local coordinates can be defined only locally and the long mode curvature perturbation enters through a mapping from local to 
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global coordinates. This coordinate transformation depends on a random field and is reminiscent of the mapping between real 
and redshift space, where the large-scale random field is given by the velocity field. It is known that redshift space distortions 
induce non-Gaussianity. 

The fact that the nonlinear GR effect is present in global coordinates, is similar to non-Gaussianities in the CMB that arise at 
recombination due solely to the intrinsic nonlinearity of relativistic perturbations Il27l429l . If the long-wavelength perturbations 
are on scales greater than the present horizon, they cannot have any physical effect, and can be rescaled away (this corresponds 
to properly redefining the background average temperature iiniini). However, we are interested in modes that are inside the 
horizon today, and we need to compare different patches of the sky modulated by a long mode ll^ \32\ . This provides an 
alternative understanding of a GR term that was missing in ll29l . compared to the expression for the squeezed CMB angular 
bispectrum obtained in ||28] |30l |33] . The angular bispectrum of CMB temperature anisotropies (and polarisation) has been 
estimated 12^ and shown to reproduce in this squeezed limit the full numerical calculations of the Einstein-Boltzmann 
system at second order IImUJTII . This intrinsic non-Gaussianity in the CMB, predicted by local rescaling arguments, could in 
principle be observed by future experiments. 

On the other hand, to compute the local abundance of halos of a fixed physical mass, we have to introduce smoothing on a 
local physical scale i? = (1 -f the corresponding local variance is independent of Cl, as argued by ifTTl . Thus the 

short-wavelength variance at a fixed physical mass scale is not modulated by the long-wavelength metric perturbation in GR. 

It is still an open question whether observable mass proxies such as luminosity, velocity dispersion. X-ray emission, SZ 
decrement, or gravitational weak lensing signal, correspond directly to the fixed physical mass defined above. Intrinsic non- 
Gaussianity in the large-scale density field in GR is real and physical. It could leave observable signatures in the large-scale 
structure of the Universe, provided we can hnd observable quantities that depend on global variables such as the distance between 
the observer and a distant patch. 
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